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Abstract. We present multi-frequency radio observational results of the quasar 
3C 48. The observations were carried out with the Very Large Array (VLA) at 
five frequencies of 0.33, 1.5, 4.8, 8.4, and 22.5 GHz, and with the Multi-Element 
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Radio Linked Interferometer Network (MERLIN) at the two frequencies of 1.6 
and 5 GHz. The source shows a one-sided jet to the north within 1", which 
then extends to the northeast and becomes diffuse. Two bright components (N2 
and N3), containing most of the flux density are present in the northern jet. 
The spectral index of the two components is ajv2 ~ —0.99 ± 0.12 and qjv3 ~ 
—0.84 ± 0.23 (S oc v a ). Our images show the presence of an extended structure 
surrounding component N2, suggestive of strong interaction between the jet and 
the interstellar medium (ISM) of the host galaxy. A steep-spectrum component, 
labelled as S, located 0.25" southwest to the flat-spectrum component which could 
be the core of 3C 48, is detected at a significance of > 15<r. Both the location and 
the steepness of the spectrum of component S suggest the presence of a counter-jet 
in 3C 48. 
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2 The radio counter-jet of the QSO 3C 48 

1. Introduction 



The radio source 3C 48 (B0134+329) is a 16th magnitude quasar, with redsift z = 0.368 
(Spinrad et al. 19851) . It was the first QSO identified at optical wavelengths (Smith et 
al. I1961f) . and the second identified one (after 3C 273) at radio wavelengths (Schmidt 
1963}. The source shows unexpectedly strong emission in the far IR bands, suggesting 
that it is located in a gas-rich environment (Chatzichristou et al. 1999). The variability 
measured with ROSAT suggests that at least part of the X-ray emission (principally the 
softest) arises from the region of the accretion-disc or from a beamed blazar component 
(Worrall et al. [2004). The image obtained with University of Hawaii's 2.2 meter optical 
telescope at Mauna Kea Observatory is dominated by a previously unsuspected high 
surface brightness region, which was referred to as "3C 48A", centered ~ 1" (1" ~ 4.1 
kpc of projected linear size x ) to the north of the QSO nucleus ( Stockton & Ridgway 
119911 Hook et al. I1994|l . This region was identified as the potential second nucleus of 
a galaxy in the process of merging with the host galaxy of 3C 48 (Zuther et al. I2003J) . 
Evidence indicates that the nuclear activity in 3C 48 is triggered by the merger process. 

3C 48 is classified as a Compact Steep Spectrum (CSS) source based on its steep 
radio spectrum and its small angular size (Fanti et al. 1985). The spectrum is rea- 
sonably straight from ~ 178 MHz up to 2.7 GHz, with a integrated radio index 
a ~ —0.83 (Peacock et al. 1198211 . It shows a low-frequency peak at about 100 MHz 
(http://nedwww.ipac.caltech.edu). The quasar 3C 48 does not show significant vari- 
ations with time and have a simple spectrum at MHz/GHz wavelengths (Baars et al. 
I1977II . and it is usually used as the primary flux density calibrator in interferometer ob- 
servations owing to its compactness and brightness. While Ott et al. 1994 found that the 
source could have become brighter from 1960s to 1990s, the increase in flux density is 
enhanced at shorter wavelengths, as is consistent with the activity in the core (Wilkinson 
et aLEHSinjl. 

Most of the radio emission occurs on a scale smaller than about 0.6" and additional 
emission extends to ~ 6" northwest to the nucleus, with a clumpy morphology (Hartas 
at al.[1983. Woan 1992). At sub-arcsecond resolution, the source exhibits an elongated 
structure in the north-south direction. VLBI observations showed a disrupted radio mor- 
phology (Wilkinson et al. 1991). The active nucleus, which is identified from its inverted 
radio spectrum (Simon et al. 1990), is weak at radio frequencies and is buried deeply 
within the host galaxy. A one-sided jet stretches out of the core towards the north ~ 1" 
within the nucleus of the host galaxy. The northern jet is highly dispersed at 0.05"- 
0.5", and its diffused morphology has been interpreted as the result of the jet outflow 
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Table 1. Observational information and image parameters 
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a the VLA observations at 0.33 and 1.5 GHz were carried out in spectral line mode, and 
recorded with narrow bandwidth; 
6 on-source integration time. 



interaction with the dense ISM in the host galaxy (Wilkinson et al. 1991). There is also 
morphological evidence for a recent merger event, such as a possible double nucleus and 
a tidal tail extending several arcseconds to the northwest (Chatzichristou et al. 119991 
Canalizo et al. lzOTUl Zuther et al. l2U03)l . 

In this paper, we present results of the multi-frequency VLA and MERLIN observa- 
tions of 3C 48. The observations and data reduction are described in the next section. 
The analysis of the images is given in Section 3. Section 4 is a discussion. A summary is 
given in the last section. 

2. Observation and data reduction 

Quasar 3C 48 was observed as a flux density calibrator with the VLA on Jan. 29, 1999 
(observational code is AH635) at 8.4 and 22.5 GHz, on Dec. 1, 2000 (AH721) at 8.4 
and 22.5 GHz, as well as on June 17, 2003 (AZ143) at 0.33, 1.5, 4.8 and 8.4 GHz. The 
observations were carried out in snapshot mode. 3C 48 was observed at 4.8, 8.4 and 22.5 
GHz with the VLA in the C configuration with a bandwidth of 50 MHz in two IFs, while 
the observations at 0.33 and 1.5 GHz were done in the spectral line mode in order to 
reject radio frequency interference (RFI). Two epochs of archival MERLIN data at 1.6 
and 5.0 GHz were also used in this paper. The observational information of the VLA 
and the MERLIN is listed in Table 1 in a time series. Column 1 is the serial number. 
Columns 2 and 3 are the epoch and the array configuration. The observing frequency, 
bandwidth, and on-source time are listed in columns 4, 5, and 6, respectively. 



4 The radio counter-jet of the QSO 3C 48 

The raw VLA data were calibrated and corrected for the phase and amplitude er- 
rors using the Astronomical Image Processing System (AIPS) of the National Radio 
Astronomy Observatory (NRAO). The MERLIN data were calibrated with the suite of 
D-programs HThomasson 1 986). 

The further process in the reduction of the VLA and MERLIN data, including editing, 
self-calibration and image analysis, was performed using both the AIPS and the DIFMAP 
packages (Shepherd, Pearson & Taylor 1994). In order to do analysis of the jet structure, 
we fitted Gaussian models to the self-calibrated uv data-sets. The task MODELFIT in 
the DIFMAP package was used for all data-sets except for the 5 GHz MERLIN, where 
we used the task JMFIT in AIPS in order to take out the extended emission prior to the 
model fitting. We found that 80-90% of the total flux density of the source is dominated 
by a few compact components. Elliptical Gaussian models were preferentially used. If 
elliptical Gaussian models failed in the least-square fitting, then a circular Gaussian was 
used. The errors in the model fitting parameters are estimated based on the method 
described by Fomalont l|1999|l . The uncertainties of flux density, originating from the 
calibration precess, are estimated within a range between 5% (for VLA) and 10% (for 
MERLIN) . The uncertainty of the component position is proportional to the FWHM of 
the synthesized beam and is inversely proportional to the ratio of peak signal to noise 
(Fomalont H5331l . 

3. Analysis 

3.1. Images 

We imaged the source using the data-sets listed in Table 1. The results show that the 
source is resolved when the resolution is higher than 0.5". Fig. 1 shows four images at 1.6, 
5, 8.4 and 22.5 GHz with sub-arcsec resolutions. The detailed structure of the source is 
shown in the images from Fig. la to Id based on an order of increasing resolution. These 
images are made using uniform weighting. Since most of the observations lack nearby 
phase calibrators, the bright, compact component N2 was used as the phase reference in 
the self-calibration and imaging process. The parameters of each image are summarized 
in Table 1. Columns 7 and 8 give the peak flux density and the off-source r.m.s. noise in 
each of the images; column 9 gives the size of the synthesized beam; column 10 indicates 
the label of the images that are present in Fig. 1. Each of the major components in the 
images are fitted in least square to a Gaussian model. The results of model-fit to the 
calibrated data are listed in Table 2. Column (1) is the serial number of the observations, 
which is the same as that in Table 1; column (2) is the observing epoch; column (3) is 
the observing frequency; column (4) is the component label; column (5) is the integral 
flux density of Gaussian component; columns (6) and (7) are the peak position of the 
Gaussian component in RA^igso and Dec^igso coordinates; columns (8), (9), and (10) 
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Fig. 1. The multi frequency total intensity images of 3C 48. The peak flux density, beam 
size, off source r.m.s level are referred to in columns 7-9 in Table 1. Contours levels for 
Figs la, lb, and lc are 1.1, 5.0, and 1.2 mJy/bx(— 1, 1, 2,..., 1024), respectively, and 
those for Fig. Id are 1.1 mJy/bx(-l, 1, 4, 16, 32, 64, 128, 256, 512). The lowest contour 
levels are generally 3 times of the r.m.s levels. 

are the sizes of major and minor axes and the position angle. Previous observations show 
that the source is dominated by steep spectrum components and the radiation from 
the core is rather weak and is difficult to identify at radio frequencies. We note that 
the phase center in each image is closer to the brightest jet component N2 (RAbi9 50 : 
01 /l 34 m 49.826 s , DEC S i 950 : +32°54'20.260") than the radio core. The core component 
and its location will be discussed in Section 4.1. 

Fig. la is the VLA image at 8.4 GHz observed at the epoch of 2003. The image 
shows an elongated structure in the NE-SW direction, which includes 4 components 
(Nl, N2, N3, and S). The main jet is dominated by the components Nl, N2, and N3. 
Component S has been detected for the first time at a significance level above 15<r, at 
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Table 2. Parameters of gaussian components 
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RA - 01 /l 34 m 49.810 s and DEC - 32°54'19.626". An independent VLA image at 8.4 
GHz, made from the observation at the epoch 2000, confirms the detection. The epoch 
2000 image at 8.4 GHz is not shown, while the parameters derived from model fitting 
are listed in Table 2. 

The two major components (N2 and N3), the northern extended component (Nl), 
and the weak southern component (S) are also detected in the 1.6 GHz MERLIN image 
(Fig. lb). A Gaussian fit failed in fitting the diffuse component Nl, visible in the MERLIN 
1.6 GHz image. Component S is confirmed at 12ct, and can be fitted with a Gaussian of 
0.28" x 0.23" (PA = 37.6°) with peak brightness of 26 mJy per beam. 
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Fig. 2. a) Simultaneous spectral index distribution (grey) between 8.4 and 22.5 GHz 
at the epoch 2000.92, superimposed on the 22.5 GHz VLA image, b) Spectral index 
distribution (grey) between 1.6 and 5 GHz, superimposed on the 5 GHz MERLIN image 
at the epoch 1992.45; Restored beam of the contour images are 0.20" x 0.20" and 0.23" x 
0.19", p.a. 45°; The contour levels are 1.3mJy/b x (-1,1,4,64,256) and 3.5mJy/b x 
(—1,1,4,64,256), respectively. The scale of the spectral index is shown in each panel. 
The darkest region corresponds to the flattest spectrum. 



Both components Nl and S appear to be completely resolved in the high resolution 
images at 22.5 GHz (Fig. lc) observed with VLA and at 5 GHz observed with MERLIN 
(Fig. Id). At the high resolutions, the component N2 and N3 are further resolved into 
subcomponents, labelled as N2a and N2b, N3a and N3b, respectively. 



3.2. Spectral index 

The multi-frequency observations allow us to study the spectral index distribution of 
3C 48 at sub arc-second resolutions. Fig. 2a shows the image of the spectral index 
distribution (grey scale) derived from the VLA images at 8.4 and 22.5 GHz at the epoch 
of 2000.92, with the contours of the radiation intensity at 8.4 GHz superimposed. The uv- 
data at 22.5 GHz is tapered to match the synthesized beam of the data at 8.4 GHz. Both 
images at 8.4 and 22.5 GHz are convolved to the same resolution and are aligned with 
respect to the common phase center. With the same method, we obtained the spectral 
index image (Fig. 2b) with the MERLIN data sets. Fig. 2b shows the spectral index image 
(grey scale) between 1.6 and 5 GHz, with the MERLIN 5 GHz contours overplotted. The 
dark end of the grey scale corresponds to the flat spectrum. Most of the source shows a 
steep spectrum. The spectrum is steepest in the northern jet (component N2), while it 
becomes flatter to the south. 

In order to image the distribution of the spectral index across the source, we plotted 
the spectral index along the jet axis at position angle of 8.5°. The slice spectral indexes 
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Fig. 3. Slice plot of the spectral index at PA 8.5°corresponding to Figs. 2a and 2b, re- 
spectively. The slice center for a) is RA 01 h 34 m 49.827 s , DEC +32°54'20.250" and for b) 
RA 01 /l 34 m 49.827 s , DEC +32°54'20.268", respectively. 

are shown in Figs. 3a and 3b, corresponding to Figs. 2a and 2b, respectively. There is a 
component centered at RA - 01 h 34 m 49.820 s ± 0.020 s , DEC - 32°54'19.740" ± 0.013", 
marked with a cross in Figs.2a&b, showing rather flat spectrum. 

The spectral index of components N2, N3 and S are estimated based on the model 
fitting results. Fig. 4 shows the power-law fitting in logarithmic coordinate. The spectrum 
of component N2 (a N2 ~ —0.99 ± 0.12) is slightly steeper than that of component 
N3 (a/v3 ~ —0.84 ± 0.23); and component S shows also a steep spectrum with a s ~ 
—0.70 ± 0.14 between 1.6 and 8.4 GHz. The results are consistent with those shown in 
Figs. 2 and 3. A mean spectral index a ~ —0.87 ± 0.10 for the overall source (3C 48) 
is estimated from the total flux densities at the frequencies between 0.33 and 22.5 GHz 
(Fig. 4). 

4. Discussion 

4.1. The position of the nucleus 

In Section 3.2 we showed that the flat-spectrum component centered at RA ~ 
01 /l 34 m 49.820 s ± 0.020 s , DEC - 32°54'19.740" ± 0.013", is most likely the nuclear core 
in 3C 48. The flat-spectrum core is likely to be buried in the strong emission at the 
centimeter and longer wavelengths. A detection of the active core was also claimed by 
Simon ct al. (1990). The core was identified at the southern end of the source in their 
22.5 GHz image. Wilkinson ct al. (1991) suggested that the radio core detected by Simon 
et al. 1990 is associated with the southernmost component in their high resolution image. 
Based on the analysis of our data, we find that the position of the flat-spectrum core in 
our image is 0.15" south of that identified by Simon et al. 1990 and Wilkinson et al. 1991. 
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Further higher resolution and higher dynamics range images at millimeter wavelengths 
arc needed to unambiguously identify the nucleus in 3C 48. 



4.2. A counter-jet in 3C 48 

We have distinctly detected an extended component S ~ 0.25" southwest of the flat- 
spectrum core, which is the putative nucleus of 3C 48. Component S shows a steep 
spectrum, i.e. a s ~ — 0.70± 0.14. Due to its location and the steep spectrum nature, we 
suggest that S is a counter-jet. 

Alternatively, a multi-particle model has been successfully developed and has been 
applied to the host galactic nucleus of 3C48 (Scharwachter et al. I2004J) . suggesting a 
counter tidal tail extending to the southwest in front of the main body of the 3C 48 host. 
The detection of component S could be a supportive evidence for the counter tidal tail 
in this quasar. 

4.3. Interaction with ISM 

CSS sources usually have sub-kpc sizes with complex radio morphologies. The nature of 
CSS sources is not completely understood yet and two scenarios are suggested. In the 
youth scenario, CSS sources are small and could reflect an early stage in the evolution 
of radio sources (Phillips & Mutel 1982, Fanti et al. I1995fl . According to the frustration 
scenario, , the unusual conditions in the interstellar medium of their host galaxies, such as 
a higher density and/or the presence of turbulence, inhibit the radio source from growing 
to larger sizes (van Breugel et al. I1984fl . The youth scenario is very much preferred 
with statistical research (Saikia et al. 119951 127M1 Murgia et al. lTM9ll2T)rTT)l . Despite this 
statistical interpretation, the frustration scenario may still apply to individual sources 
which can be interpreted in terms of confinement and interaction between the jets and 
ISM in the host galaxies. Numerical simulations have also shown that a dense ISM may 
substantially confine CSS sources (De Young 1992J)- The radio structure of 3C 48 has 
shown the evidence for strong interaction between the radio jet and the ISM in the host 
galaxy. Here is a working scenario. A supersonic collimated jet stretches out from the core 
to the north where a shock is produced in the interaction between jet flow and the dense 
medium. The interaction converts the kinetic energy carried by the jet to internal energy 
or turbulence which could enhance the acceleration of particles and then the radiation 
field, causing the jet to brighten. The observed north jet brightening (hot spot N2 in 
Fig.lc and Fig. Id) after losing collimation appears to be consistent with this picture. 
Based on the study of the kinematics and physical conditions of the extended emission 
line gas and their relations to the nuclear star formation, Chatzichristou (2001) found that 
the coupling between radio jet and gas is the kinematic signature of a recent interaction. 
Both Canalizo & Stockton (2000) and Chatzichristou et al. (1999) suggested that the 
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Fig. 4. Total spectrum of 3C 48 (solid line) and of the various components in the source. 

unusually high stellar velocity dispersion as well as very young stellar populations are 
possibly related to the interaction of the nuclear radio jets with the ambient media. The 
NIR (near infrared) images and spectra also suggest that reddening by several magnitudes 
in the nuclear emission can be explained by the interactions of the 3C 48 radio jet with 
the circum-nuclear medium (Zuther et al. I200 3j . The ratio of the recombination lines 
P a j3 and P a j determined from 3C 48 (Zuther et al. 2003) is about 1.6 (±0.5), which is 
consistent with the value typical for the ISM of active galactic nuclei being excited by 
the central-engine continuum emission. 



5. Summary 

In this paper, we carried out radio multi-frequency observations of 3C 48. The high 
dynamic range images exhibit the detailed structure of the source on the sub-arcsec scale. 
The position of the nucleus with a flat spectrum is determined based on the spectral index 
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distribution images. The jet components show steep spectra. We confirm the existence of 
a counter-jet in 3C 48, which shows the evidence for the first time that 3C 48 actually 
has a double-sided jet. The northern jet appears to be disrupted at 0.5"northeast to 
the core. Our analysis suggests that the initially collimated jet can be disrupted by the 
interaction between jet flow and the circum-nuclear medium. 
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